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Abstract —Time-varying renewable energy generation can re¬ 
sult in serious under-/over-voltage conditions in future dis¬ 
tribution grids. Augmenting conventional utility-owned voltage 
regulating equipment with the reactive power capabilities of 
distributed generation units is a viable solution. Local control 
options attaining global voltage regulation optimality at fast 
convergence rates is the goal here. In this context, novel reactive 
power control rules are analyzed under a unifying linearized grid 
model. For single-phase grids, our proximal gradient scheme 
has computational complexity comparable to that of the rule 
suggested by the IEEE 1547.8 standard, but it enjoys well- 
characterized convergence guarantees. Adding memory to the 
scheme results in accelerated convergence. For three-phase grids, 
it is shown that reactive injections have a counter-intuitive effect 
on bus voltage magnitudes across phases. Nevertheless, when 
our control scheme is applied to unbalanced conditions, it is 
shown to reach an equilibrium point. Yet this point may not 
correspond to the minimizer of a voltage regulation problem. 
Numerical tests using the IEEE 13-bus, the IEEE 123-bus, 
and a Southern California Edison 47-bus feeder with Increased 
renewable penetration verify the convergence properties of the 
schemes and their resiliency to grid topology reconfigurations. 

Index Terms —Linear distribution flow model, accelerated 
proximal gradient, three-phase distribution grids, PV Inverters. 

I. Introduction 

Voltage regulation, that is the task of maintaining bus 
voltage magnitudes within desirable levels, is critically chal¬ 
lenged in modern distribution grids. The penetration of re¬ 
newables, demand-response programs, and electric vehicles 
lead to time-varying active power injections and frequently 
reversing power flows. Utility-owned equipment convention¬ 
ally employed for voltage regulation, such as tap-changing 
transformers and shunt capacitors, cannot react promptly and 
efficiently enough m, ca. Hence, avoiding under-Zover- 
voltage conditions requires advanced reactive power manage¬ 
ment solutions. To that end, the power electronics of PV 
inverters and storage devices offer a decentralized and fast- 
responding alternative ll 28 l . ||30l- 

A grid operator can engage the reactive power capabili¬ 
ties of distributed generation (DG) units to minimize power 
losses while satisfying voltage regulation constraints. Being 
an instance of the optimal power flow (OPF) problem, re¬ 
active power support can be solved using convex relaxation 

Work in this paper was supported by NSF grants 1423316, 1442686, 
1508993, and 1509040. V. Kekatos is with the ECE Dept., Virginia Tech, 
Perry Str, Blacksburg, VA 24060, USA. L. Zhang and G. B. Giannakis 
are with the Digital Technology Center and the ECE Dept., University of 
Minnesota, Minneapolis, MN 55455, USA. R. Baldick is with the ECE 
Dept., University of Texas at Austin, TX 78712. Emails: kekatos@vt.edu, 
{zhan3523,georgios}@umn.edu, baldick@ece.utexas.edu 


techniques Eoi, ini. Among other centralized approaches, 
inverter VAR control is solved using convex relaxation in ifOll . 
while a scheme relying on successive convex approximation 
is devised in m. Distributed algorithms requiring commu¬ 
nication across nodes have been proposed too. A distributed 
method based on convex relaxation has been developed in ll^ . 
Upon modeling power losses as a quadratic function of reac¬ 
tive power injections, lH pursues a consensus-type algorithm. 
Control rules based on approximate models are presented 
in ll^ . and ^ developes a multi-agent scheme. Building on 
the radial structure of distribution grids, algorithms requiring 
communication only between adjacent nodes have been devel¬ 
oped based on the alternating-direction method of multipliers 

(ADMM) m, ca, El. 

To cater the scalability of DG units and the potential lack 
of communication infrastructure, local plug-and-play schemes 
are highly desirable. Given that voltage magnitudes depend on 
grid-wide reactive injections, guaranteeing voltage regulation 
constraints may be hard to accomplish via purely localized 
algorithms. In such setups, reactive power management is 
usually relaxed to penalizing voltage magnitude deviations 
from a desired value and neglecting power losses. Reactive 
power injections are adjusted proportionally to the local volt¬ 
age violations in lIZTl : see also ll^ for sufficient conditions 
guaranteeing its convergence. A similar local control strategy 
has been shown to minimize a modihed voltage regulation cost 
in QSl, while ifTrll proposes a subgradient-based algorithm. A 
control rule adjusting the inverter voltage output according to 
the reactive power flow is reported in |[T^ . The local control 
rules proposed in HTsj maintain voltage magnitudes within 
the desired range under the presumption of unlimited reactive 
power support. 

Most existing works build on a simplified single-phase grid 
model. Due to untransposed distribution lines and unbalanced 
loads though, the equivalent single-phase distribution network 
may not exist. Semidefinite programming and ADMM-based 
schemes have been applied in multiphase radial networks for 
power flow optimization ini, ini. Nonetheless, no work 
discloses how inter-phase coupling affects bus voltage mag¬ 
nitudes across phases, or how local control schemes behave in 
unbalanced grids. 

This work considers local reactive power control rules for 
voltage regulation in single- and multi-phase distribution grids. 
Our contribution is on four fronts. First, we provide a unified 
matrix-vector notation for approximate yet quite accurate 
multi-phase grid models (Section [Till. Second, Section HII] 
extends the work of mi. In mi, we developed a reactive 
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power control rule based on a proximal gradient scheme, 
and engineered modifications with superior convergence prop¬ 
erties. Here, the options of unlimited reactive support and 
diagonal scaling are considered too. In particular, numerical 
results indicate that the convergence rates attained by using 
different step sizes across buses are still significantly lower 
than those achieved by our accelerated scheme. Third, using 
a linear approximation for unbalanced multi-phase grids, we 
reveal an interesting inter-phase coupling pattern across buses 
(Section HVl) . Recall that in single-phase grids, increasing the 
reactive power injection at any bus raises the nodal voltage 
magnitudes throughout the grid. In multi-phase grids on the 
contrary, injecting more reactive power into a bus of one 
phase could result in decreasing voltage magnitudes for the 
preceding in the positive-sequence ordering phase. It is finally 
shown that in unbalanced scenarios our reactive power control 
rule converges to an equilibrium point; yet this point does not 
necessarily correspond to the minimizer of a voltage regulation 
problem. Numerical tests on distribution feeders corroborate 
the convergence properties of the novel schemes, as well as 
their resiliency to topology reconfigurations (Section fVTl. 

Regarding notation, lower- (upper-) case boldface letters 
denote column vectors (matrices), with the exception of line 
power flow vectors (P, Q). Calligraphic symbols are reserved 
for sets. Symbol ^ stands for transposition. Vectors 0, 1, 
and e„, are the all-zeros, all-ones, and the n-th canonical 
vectors, respectively. Symbol ||x ||2 denotes the ^ 2 -norm of 
vector X, while Ai(X) stands for the i-th largest eigenvalue 
of X. Operator diag(x) defines a diagonal matrix having x 
on its diagonal, whereas diag(X) is the vector formed by the 
main diagonal of X. Operator bdiag({X„}) defines a block 
diagonal matrix with X„’s as blocks. Operators Re(z) and 
Im(z) return the real and imaginary part of complex number z. 
A matrix with non-negative entries is denoted by X > 0, while 
X A 0 indicates a symmetric positive semidefinite matrix. 

II. Radial Distribution Grid Modeling 

Approximate models for single- and multi-phase grids are 
presented in this section. 

A. Modeling of Single-Phase Grids 

Distribution grids are typically operated as radial. A single¬ 
phase grid with N -\-l buses can be modeled by a tree graph 
T — {ffo, C) whose nodes Afo ■= {0, ..., N} correspond to 
buses, and whose edges £ correspond to distribution lines with 
cardinality \C\ = N. The feeder bus is indexed by n = 0, 
whereas every non-feeder bus n G Af = {1,...,A(} has a 
unique parent bus denoted by 7r„. Without loss of generality, 
nodes can be numbered such that 7r„ < n for all n G Af. For 
every bus n G Afo, let be its squared voltage magnitude, 
and Sn = Pn + jqn its complex power injection. 

The distribution line connecting bus n with its parent 7r„ is 
indexed by n. For every line n G C, let z„ = Vn-GjXn, in, and 
Sn = Pn A jQn be the line impedance, the squared current 
magnitude, and the complex power flow sent from the sending 
bus TTn, respectively. If is the set of children buses for bus 


n, the grid can be modeled by the branch flow model 0, a 

Sn — ^ ^ Sk Sn “b in^n (la) 

k&Cn 

Vn = - 2Re[z^Sn] + in\Zn\^ (lb) 

|5'„P=U7r„C (Ic) 

for all n G Af, and the initial condition sg = SfcGCo 

For notational brevity, collect all nodal quantities related 
to non-feeder buses in vectors p [pi ■ ■ ■ pat]^, q ;= 
[<?! and v := [ui ••• wat]^. Similarly for lines, 

introduce vectors r := [ri ••• tat]^, x := [xi ••• xnY, 
P := [Pi ••• Pat]^, and Q := [Qi ■■■ Qn]^- Define 
further the complex vectors s := p -b jq, z := r -b jx, and 
S := P -b jQ. According to the approximate LinDistFlow 
model, the grid is described by the linear equations 0, 0 

s = A^S (2a) 

Av = 2Re[Z*S]-agug (2b) 

where ug is the squared voltage magnitude at the feeder; matrix 
Z is defined as Z := diag(z); and A is the reduced branch-bus 
incidence matrix enjoying the following properties. 

Proposition 1 filin'). The negative of the reduced branch-bus 
incidence matrix —A and its inverse F := —A~^ satisfy: 
(pi) they are both lower triangular with unit eigenvalues; 
(p2) F > 0; and Fag = Iat 

where A = [ag A] is the full branch-bus incidence matrix. 

Equation can be now expressed as S = —F^s. 
Substituting the latter into (l2bl) . and premultiplying by —F 
yields v = 2Re [FZ*F^s] -b 'UgFag. Proposition [T] and the 
properties of the real part operator provide 0, ifTSll 

V = Rp-b Xq-b uglAT (3) 

where R := 2Fdiag(r)F^ and X := 2Fdiag(x)F^; see 
also 0 for a linear approximation model relating complex 
voltages to complex injections. It is well understood that in 
transmission grids and under regular load conditions and high 
reactance-to-resistance ratios, the nodal voltage magnitudes 
are approximately independent of active power injections. 
On the contrary, the approximate model of Q confirms that 
voltage magnitudes in distribution grids depend significantly 
not only on reactive but active injections too. The dependence 
is roughly linear with the following properties. 

Remark 1. Although F is lower triangular, matrices R and 
X are generally full. Hence, local injection deviations affect 
voltage magnitudes globally. 

Remark 2. Assuming r and x to be strictly positive, R and 
X are symmetric and strictly positive-definite by definition; 
cf. im for a more elaborate proof. Further, their inverses 
are expressed as R“^ = ^A^ diag~^(r)A and X“^ = 
iA^ diag“^(x)A. 

Remark 3. Since F > 0, it follows readily that R > 0 and 
X > 0. Hence, injecting more active or reactive power at any 
bus raises the voltage magnitudes in the entire grid. 


3 


B. Modeling of Multi-Phase Grids 


where the argument inside the real operator simplifies as 


The focus shifts next to modeling multi-phase grids. For 
ease of exposition, it is first assumed that all buses are served 
by all three phases. For this reason, system variables are now 
represented by 3-dimensional vectors. Slightly abusing the 
notation used in Section Ill-AI the complex voltages and the 
power injections at all phases of bus n here are denoted by v„ 
and s„, respectively. Similarly, the complex currents and the 
complex power flows on all phases of line n are represented by 
i„ and S„, respectively. The coupling across phases provides 
the multivariate version of Ohm’s law; 


(4) 

where Z„ = is the phase impedance matrix for line n. If 
S„ is the flow on line n seen from bus 7r„, the flow received 
at bus n is 


diag(v„)i* = S„ - diag(Z„i„)i*. 

The multi-phase generalization of ( flai l reads 

s™ = ^ Sfc - S„ -f diag(Z„i„)i* (5) 

keCn 

for all n G Af. As advocated in El, to obtain the multi-phase 
equivalent of the voltage drop equation (fTbl i. multiply (l4|i by 
the complex Hermitian of each side, and maintain only the 
diagonal of the resultant matrix; 


diag (v„v^) = diag - 2 Re 

-f diag (Zninin'^n) ■ 


diag 

( 6 ) 


The full AC model for this multiphase grid is completed by 
replicating the definition for flows ( fT^ on a per phase basis. 
Similar to single-phase grids, the model involves computation¬ 
ally inconvenient quadratic equations, but convex relaxations 
render the model tractable under appropriate conditions iflTl . 

Alternatively, one may resort to the simpler approximate 
model of ifTTl . As for single-phase grids, because Z„’s have 
relatively small entries, the last summands in the right-hand 
sides (RHS) of dill and (01 can be dropped. Regarding the 
second summand in the RHS of (0, let us further assume that 
phase voltages are approximately balanced. By surrogating v„ 
by Vna, where a. :=\1 a a = e~^~, the complex 

current vector can be roughly expressed as 


diag (Z„ diag(S* )q;)^^ = diag(Q;)Z* diag(S„)a* 

= diag(Q;)Z^ diag(a*)S„. 


The equalities follow from the properties of the diag operator; 
diag(xy^) = diag(x)y* and diag(x)y = diag(y)x. The 
approximate multi-phase model reads for a\\ n G Af 



v.„ -v„ = 2Re[z;S„ 

Z„ ;= diag(a*)Z„ diag(Q;). 


(8a) 

(8b) 

(8c) 


Building on the grid model of El, we express (01 in a 
matrix-vector form and study the involved matrices to better 
understand voltage regulation schemes. Heed that system 
variables can be arranged either per bus or per phase. For 
example, the squared voltage magnitudes can be stacked as 



Va 


Vl 

V := 

Vf, 

or V := 



Vc 




where ;= [vf ... vf^]^ for f G {a,b,c}. Likewise, 
injections can be arranged in s or s, and flows in S or S. 
It can be easily verified that the aforementioned vector pairs 
are related by 

V = Tv, s = Ts, S = TS (10) 


for a common permutation matrix T compactly expressed as 


T ;= 


I3 ® 

I3 ® ejf 


( 11 ) 


where e„ is the n-th column of Ijv. Being a permutation 
matrix, T satisfies = T^. Algebraic manipulations 

postponed for the Appendix show that voltage magnitudes in 
multi-phase grids are related to nodal injections as follows; 

Proposition 2. Based on (01, it holds that 


V = Rp + Xq + wol3W (12) 


i* Ri diag(S„)a* (7) 

and the outer product can be thus replaced by 

diag(S„). Let v„ be the vector of per-phase squared 
voltage magnitudes for all three phases on bus n 


where the involved matrices are defined as 

R ;= 2Mbdiag({Re[Z„]})M^ (13a) 

X ;= 2Mbdiag({Im[Z„]})M^, and (13b) 

M ;= T(l 3 0 F)T^ (13c) 


= diag(v„v„ ). 


Then, equation 0 permits the approximation 

v,r„ - v„ = 2 Re [diag [olol^ diag(S„)Z; 


and matrices Z„ have been defined in (I8cl i. 

According to (ISal l. power injections and flows are approxi¬ 
mately decoupled across phases. Nonetheless, Proposition 0 
asserts that squared voltage magnitudes depend on power 
injections from all phases. Building on the approximate models 
of this section, voltage regulation schemes are developed next. 
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III. Schemes eor Single-Phase Grids 

Posing voltage regulation as an optimal power flow instance 
leads to a constrained optimization problem. Given that voltage 
regulation constraints couple reactive injections across the 
grid, developing localized solutions becomes challenging. To 
derive such solutions, the voltage regulation goal is relaxed 
here and posed as the generic minimization problem 

min/(q)+c(q) (14) 

q6Q 

where /(q) is the cost of squared voltage magnitudes v devi¬ 
ating from their nominal value uqI; c(q) models the potential 
cost for reactive power compensation; and Q is the feasible set 
of reactive injections. Since reactive power injections by DG 
inverters can be adjusted in real-time, whereas utility-owned 
voltage regulating equipment responds typically at a slower 
pace (e.g., every few minutes or hourly), the latter choice is 
assumed hxed to a value and it will not be considered here. 
Particular instances of the generic setup in (fl4l i are instantiated 
next for single-phase grids. 

A. Unconstrained Reactive Support 

A viable voltage deviation cost is /i(q) := i||v — uollH- 
This cost function tries to keep squared voltage magnitudes 
close to the nominal value vq. Assuming c(q) = 0 and Q = 
problem (fT4l i boils down to the unconstrained quadratic 
program 

mini||Rp + Xq||2. (15) 

q 

Localized voltage regulation schemes assuming unlimited 
reactive power support have also been considered in ll3?l 
and ll2^ . Obviously, since X is invertible, problem (flsT l has 
the unique minimizer q* = — X“^Rp that achieves perfect 
voltage regulation v(q*) = uqI. Finding q* requires knowing 
the active injections p over all buses. Using the structure of 
R and X, vector q* can be alternatively expressed as 

q*=ATdiagQ^|)p. (16) 

The entry q* is a linear combination of the active powers 
flowing in and out of bus n with the related r„/a;„ ratios as 
coefficients. Although denotes the power flow seen from 
the sending end of line n, the receiving end will measure 
approximately — due to the small loss assumption. The 
minimizer (Ell can be found in a localized way only if bus n 
measures power flows on incident lines. 

Alternatively, a gradient descent scheme would iteratively 
update reactive injections over time t as 

q‘+i = q* - /ig* (17) 

where p, > 0 is a step size, and g‘ = X^(Rp + Xq*) = 
X(v* — uqI) is the gradient of /i(q) at q*. Unfortunately, 
such a scheme cannot be implemented in a localized fashion. 
However, the next proposition proved in the Appendix asserts 
that the rule 

q‘+i = q* - p(v* - uol) (18) 


Proposition 3. If p € (O, 2Aniin(X)/A^3^x(X)), the rule of 
(El converges to the minimizer of (E)- 

The descent rule in (El l scaled by a diagonal matrix has 
been shown to converge for a more detailed model in ll^ . 

B. Constrained Reactive Support 

Solving (Ell may be of little practical use: Distributed 
generation units may not be installed on every bus and their 
reactive power resources are limited. Moreover, the power 
electronics found on a PV at bus n have finite apparent power 
capability s„, which limits (?„ depending on the current active 
injection (solar generation) as < s^. In reality, q is 

constrained to lie within the time-varying but known box 

Q := {q : |g„| < := s/si-pi > 0 Vn}. 

Buses with reactive power support can be obviously modeled 
by selecting their associated limits as = 0. 

In this practically pertinent setup where the voltage regula¬ 
tion problem in (fl4l l is constrained as q G Q, one could try 
implementing the projected version of (fEll . that is 

q*+i = [q*-^(V-uol)]g (19) 

where [x]q := argmin^gg ||x — z ||2 denotes the projection 
operator on Q. Unfortunately, this seemingly meaningful con¬ 
trol rule is not guaranteed to converge Q. 

A localized voltage regulation scheme can be obtained via 
a different voltage deviation cost /(q). As advocated in ifTSll . 
a meaningful choice is the cost 

/ 2 (q) = i||Rp + Xq||^_, (20) 

with the rotated norm dehned as Ijxljg := x^Bx for 
B 0. As proved in 121], the cost in (l20t is equivalent to 
/ 2 (q) = i En=i ■ Although minimizing / 2 (q) over 

q G Q penalizes scaled voltage magnitude deviations between 
adjacent buses, obviously, it does not guarantee that voltages 
will lie within any desired range. Nevertheless, / 2 (q) has the 
important feature that its gradient 

V/ 2 (q) =Rp + Xq = v(q)-uol (21) 

equals the deviation of squared voltage magnitudes from the 
nominal, and its n-th entry can be measured locally at bus n. 

To deter engaging PV power inverters for negligible voltage 
deviations, a reactive power compensation cost c(q) should 
be also considered. Given that negative and positive reactive 
power injections are equally important, a reasonable option for 
voltage regulation would be solving the problem 

min /i 2 (q) = / 2 (q) + C 2 (q) (22) 

qeQ 

where C 2 (q) := J2n=i^n\qn\- Again, due to the strong 
convexity of / 2 (q), problem (l22li has a unique minimizer in 
Q. As shown in IImI . the minimizer of (l22l l can be found via 
simple proximal gradient iterations: At iterate t, each bus n 
measures the quantity 

J/n := - ^'o) 


converges to q* for an appropriately small p. 


(23) 
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Qn,t 



for a step size fj, > 0. Voltage magnitude deviations — vq 
are assumed to be acquired without noise. It then updates its 
reactive power injection according to the control rule 




The operator S^{y,c) is defined as 


(24) 




+q 

y- Hc 
0 

y + yc 


-q 


, y>q + y,c 
, fic < y <q + yc 
, — yc < y < yc 
, -q- yc<y < -yc 
, y < -q- yc 


(25) 


and is shown in Figure[T] Apparently, in the absence of reactive 
power compensation cost, that is when c — Q, the operator 
5^(j/,0) simply projects y onto [—q,q]. 

If the step size is selected as y G (0,2A“^x(^))’ the 
control rule of (l2Tll- (l24l l converges to the minimizer of (|2^ 
ll25l . Actually, with the optimal step size y = A“t^(X), the 
convergence rate is linear, but proportional to the condition 
number At(X) := Aniax(X)/Aniin(X) of matrix X. It is worth 
stressing that k(X) can be relatively large: the Southern 
California Edison 47-bus grid and the IEEE 34-bus benchmark 
exhibit k(X) = 1.6x 10^ and 5.5x 10"*, respectively ifTSl , BTl . 
Ways to improve the convergence speed are reviewed next. 


C. Accelerating Voltage Regulation Schemes 

A way to speed up proximal gradient schemes is via 
diagonal scaling. In detail, the reactive injection variables q 
can be transformed to q as q = D^/^q for a diagonal matrix 
D 0. Instead of (l22l l. we can solve the equivalent problem 

min hzfD^/^q) 

through the proximal gradient method. The update for the 
transformed variable q‘+^ is found as the solution to 

min C2(Di/2q) -f ^\\q- {q* - yU^/'^Vf2i'D^^‘^q*))\\l 

and the convergence rate now depends on k 

Translating the iterations back to the original variables q yields 

q*+i:=argminc2(q)-|-^||q-(q‘-/iDV/2(q*))||D-i (26) 


which is separable across buses as 

1 2 

= arg min d„c„|q„|-f — - uM (27) 

q^<qn<qn Zy 

where := — ydn{vl^ — vq). It can be easily verified that 

the minimizer of (|26]) is provided by the control rule 

ql^^ ■= 'SldJul.^dnCn) (28) 

where dn is the n-th diagonal entry of the scaling matrix D. 
The control rule in (l28l l corresponds to the rule of (l2^ - (l24l) 
with the step size y adjusted to ydn per bus n. The scaling 
matrix D could be selected to minimize n A 

relatively simple choice for D is to assign vector diag(X) on 
the main diagonal of D~^. This option sets the diagonal entries 
of D^/^XD^/^ to unity. To meet faster voltage regulation 
rates, the accelerated voltage regulation scheme of ED is 
reviewed and simplified next. 

In ED, we derived an accelerated proximal gradient scheme 
by adapting Nesterov’s method EH. The corresponding con¬ 
trol rule was shown to be a simple modification of EH (see 
IITTI for details): 

(29) 

where variable is updated using the from (l23t as 

yl:={l + f3t)yl-f3tyl^ (30) 

with fit = for all t > 1. Compared to (I241 l. the 
control rule in EU introduces memory in calculating y!f as 
a linear combination of y!f and The linear combination 
coefficients depend on the time-varying fit that converges to 
1. Note that the formula for fit has been simplified from the 
one used in ETl Eqs. (20)-(21)]. 

If the step size is selected as y = A“ax(^) ^od the sequence 
fit is reset to zero every 2 y^k(X) iterations, the reactive 
control rule of (l29l l converges linearly to an e-optimal cost 
value within iterations. Eor grids with high 

«:(X), this modihed scheme offers accelerated convergence. 
Numerical tests show that even without resetting the sequences 
and without knowing precisely Aniax(X), the rule in EH) offers 
superior performance over both rules EH and (l28l l. 

Remark 4. The IEEE 1547.8 standard suggests the following 
reactive power injection rule for the DG at bus n im 

q^^'^ = Sn{vo - vlf) (31) 

where the function Sn is an increasing piecewise linear func¬ 
tion similar to the one shown in Eigure [D Comparing OTI) 
with the control rules of EH or (l28T l. suggests that all control 
rules have similar computational requirements. Nevertheless, 
the rule in ED can exhibit oscillations as reported in lfT4l . 
while the schemes presented here exhibit well-understood 
convergence properties. 

IV. Schemes for Multi-Phase Grids 

Distributed generation and demand-response programs can 
lead to unbalanced conditions. This section aims at general¬ 
izing the schemes of Section HII] for multi-phase grids. The 
problem of voltage regulation in multi-phase grids can be 
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posed as in (fT4l t: yet now the dependence of squared voltage 
magnitudes on nodal injections is governed by the model 
in dnii. Before devising voltage regulation schemes, critical 
properties of the involved matrices are studied first. 

Let us focus on the 3x3 complex matrices Z„ defined in 
(| 8 ]l. Let + jx’^i be the (i,j)-th entry of Z„. From 

the symmetry of Z„ and the identity = a* = j ^7 

matrix Z„ becomes 


Z„ = diag(a*)Z„ diag(Q:) = 


vii 


az, 


21 


a z, 
^22 


12 


az. 


13 


a z: 


31 


az: 


32 


a z. 

.,33 


23 


Therefore, matrix Im[Z„] can be decomposed as 


Im[Z„] 

= X„ 

+ Rn 




1 

o 



•''n 

^13 

x„ 


^,12 


^23 


z 


^13 

™23 

•''n 

2x33 

Rn 


0 

j,12 

n 

_„12 
' n 

0 

„13 
' n 

_„23 
' n 


z 


1 

„23 
' n 

0 


(32a) 

(32b) 

(32c) 


where X„ is a symmetric matrix (X„ = X^) associated to 
the reactive part of Z„, and R„ is an anti-symmetric matrix 
(R„ = —R^) depending on the resistive part of Z„. Using 
(l32T i. the next fact holds: 


Proposition 4. If X„ is strictly diagonally dominant with 
positive diagonal entries, then Im[Z„] 0. 

To prove Proposition @1 it suffices to show that the symmet¬ 
ric component of Im[Z„] is strictly positive-definite. If mutual 
and self-reactances satisfy 2 x“ > l L '^ben 

X„ is diagonally dominant with positive diagonal entries, and 
thus, X„ >- 0. Due to the structure of distribution lines, the 
aforementioned inequalities are satisfied in general. 

The decomposition in (l32l l carries over to X in (I13bb as: 


X = X, + X, (33) 

where X^; := 2Mbdiag({X„})M^ is symmetric, and X^ := 
2M bdiag({R„})M^ is anti-symmetric. Because M is in¬ 
vertible, matrix X is positive-definite if and only if matrix 
bdiag({X„}) is, hence leading to the corollary: 

Corollary 1. If Im[Z„] is strictly diagonally dominant with 
positive diagonal entries for all n, then X ^ 0. 


A. Inter-Phase Coupling 

We next elaborate on how bus voltage magnitudes are 
affected by reactive power injections. According to (fT^ . vector 
V is linearly related to reactive injections q via the matrix 

X := T^XT (34) 

= 2(13 ® F)T^ bdiag({Im[Z„]})T(l 3 0 F^). 

The effect of reactive power injection to the squared voltage 
magnitude is described by the {i,j) entry of X. Let entry 


i correspond to phase fi of bus Ui, and entry j to phase 
of bus rij. It can be verified that 

^-i,j ~ ^ 0 ) 

= 0 f„J^T^ bdiag({Im[Z„]})T(e 0 ^ 0 f„J 

where f^ is the fc-th row of matrix F. By the definition 
of T in (fTTT i. the products T(e 0 . 0 f„J and T{e^^ 0 f„j) 
can be expressed as f„^ 0 e^. and f„^. 0 respectively. 
Exploiting the structure of bdiag({Im[Z„]}) and since F is 
lower triangular, the entry X^^ can be expressed as 

min{ni ,nj } 

^i,j = (35) 

fc=i 


Recall that F has non-negative entries, while for overhead 
distribution lines the parameters and appearing 

in (I32bl l- (l32cl l are typically non-negative. According to dTSl l. 
three cases can be distinguished: 

(cl) If = fj, X,,, = > 0. 

Thus, as in single-phase grids, injecting more reactive power 
into a bus raises voltages at all buses in the same phase. 

(c2) When G {(a, &), ( 6 , c), (c, a)}, it follows that 




min{ni,nj} ( 6 ^ 


Thus, injecting more reactive power into a bus decreases the 
voltage magnitudes at all buses of the preceding phase in the 
positive sequence ordering. 

(c3) If G {(a,c), ( 6 ,a), (c, &)}, then X 

_ 1 ^min{ni.n,} / ^0^ _ /o 0^0) p p , 

2 2.^fc=l y^k vO'k j ^ Ui.kl^ nj,k- 

effect of reactive power injections into one phase to the voltage 
magnitudes of the following phase depends on the differences 


< 0 . 


Thus, the 


appearing in the last sum. Actually, if every 


bus serves all phases and since 




- \/3r 


< ^k 


ipitpi 


+ 


the effect of one phase to the following phase is less 
significant than its effect on the previous phase. 

An illustration of the influence patterns across phases is 
shown in Figure |2] Evidenced by the previous analysis and 
Fig.El decreasing reactive injections in phase b to cater over¬ 
voltage conditions on phase b would aggravate possible over¬ 
voltage problems on phase a. In this context, voltage regulation 
becomes even more challenging in multi-phase grids. 


B. Unconstrained Reactive Support 

Similar to single-phase grids, let us first consider the 
simplest voltage regulation scenario: The voltage deviation 
cost is /i(q) = ^||v — wolaArlli; there is no reactive power 
compensation cost c(q), and reactive power is unconstrained: 

mini||Rp-f Xq|| 2 . (36) 

q 

Again, the positive-definiteness of X [cf. Corollary. [T] guar¬ 
antees the uniqueness of the minimize!' in (l36l l. Contrary to 
the single-phase case though, the minimize!' of (l36l l cannot be 
found locally even if bus n measures all phase flows on its 
incident lines [cf. (ITSI ll because of the structure of R and X. 

Similar to single-phase grids, a gradient descent scheme 
cannot be implemented locally. However, the simple update 

q‘+i = q‘-/i(v‘- !;ol) 


(37) 
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Fig. 2. Matrix X (normalized to unity maximum entry) relating voltage 
magnitudes v to reactive power injections q for the IEEE 13-bus grid 
benchmark depicts the coupling across phases. 


still constitutes a convergent rule: 

Proposition 5. If G A^^rTx^x) )’ (l37l i 

converges to the minimizer of (I361 l. 

The key point here is that albeit X is not symmetric, it is 
still positive-definite under the assumptions of Corollary [1] 
Hence, rule (|J7] | remains a valid descent direction for the 
unconstrained problem in (|36] |. 

C. Constrained Reactive Support 

Reactive power injections are typically constrained in a 
feasible set Q, and there may also be some reactive power 
compensation cost C 2 (q). In that case, the update in (iJTl l is not 
practical. Recall that for single-phase grids with constrained 
reactive resources, the original voltage magnitude deviation 
function /i(q) = i||v — uol||i was replaced by a rotated 
Euclidean norm of the voltage deviations, namely / 2 (q) = 
i||X“^/^(v — uol)|||, which resulted in localized updates. 
The choice of /i(q) would fail yielding localized solutions 
in multi-phase grids too. Although X in multi-phase grids is 
positive-definite, the lack of symmetry does not allow us to 
extend the approach with / 2 (q)- 

Nonetheless, let us study the behavior of the control rule of 
(l24t under unbalanced conditions. Assume that the DG unit 
at each bus performs the control rule of (l23]) - (l24l) that can be 
equivalently expressed as (cf. ED) 

y‘= q* - p(v‘- uol) (38a) 

q‘+i = prox^^^_Q[y*] (38b) 

where the proximal operator is defined as 

[y] — arg min pc 2 (w) -h i||w - y\\l (39) 

wt y 

Compared to the single-phase grid case, the major difference 
is that now v is related to q according to the model in (fTTl i. 


The iterates produced by ( l38T l satisfy: 

- q*ll2 = II proxies.eiri] - pi'ox^c2.e[y‘"^]ll2 
<l|y‘-y‘-'ll 2 
= 11 (I-pX) (q‘-q*-l) ||2 
<||I-pX||2||q‘-q*-i||2. 

where the first inequality follows from the non-expansive 
property of the proximal operator (cf. Il25l Prop. 5.1.8]); the 
equality from (I38al) : and the last inequality from the definition 
of the maximum singular value. If p is selected such that 
||I—pX ||2 < 1, then (I38I) constitutes a non-expansive mapping 
and it therefore converges to the equilibrium point defined by 

q* = prox^^^ glq* - At(v(q*) - VqI)] (40) 
or, from the definition of the proximal operator in ( |39] |. by 
q* = arg min /iC 2 (w) + i||w - [q* - p.{v{q*) - Uol)]||i 

wGQ 

Vector q* is thus defined as the minimizer of an optimization 
problem, and it cannot be expressed in closed form. Of course, 
the equilibrium point q* does not necessarily coincide with the 
minimizer of any voltage regulation optimization problem. The 
next step size range guarantees ||I — /rX ||2 < 1, and therefore 
convergence of (l3Ft (see the appendix for a proof): 

Proposition 6. Let UAU^ be the eigenvalue decomposition 
o/XX^. If p G (0,Amin(A”^/2UT(X + XT)UA”i/2))^ 
then III — /rX ||2 < 1. 

Practical distribution grids do not have all phases at all 
buses. The previous analysis carries over to such cases, if 
the related entries of Z„’s and the corresponding (re)active 
power injections are set to zero. For the eigendecompositions 
of XX^ and X -|- X^, their rows and columns related to 
non-existing bus-phase pairs are simply removed. 

V. Numerical Tests 

The voltage regulation schemes presented earlier are evalu¬ 
ated using the IEEE 13-bus feeder, the IEEE 123-bus feeder, 
and a Southern California Edison (SCE) 47-bus feeder M, 
ED- Solar generation data from the Smart* project and from 
August 24, 2011 are used ||5l. Unless otherwise stated, active 
power injections are kept fixed over the reactive control period, 
and PV reactive injections are initialized to zero. 

A. Single-Phase Grids 

Single-phase grids are obtained upon modifying multi-phase 
grids as in ED- The global minimizer of (fTTl i is obtained using 
MATLAB and it serves as a benchmark. The first experiment 
simulates an over-voltage scenario on the IEEE 13-bus grid 
depicted in Fig.|2 The IEEE 13-bus grid exhibits «:(X) = 716. 
Simulating a severe over-voltage violation in the transmission 
network, the feeder voltage magnitude is set to 1.07 p.u., while 
the voltage regulator is removed from the system. A 52% PV 
penetration level is assumed for all buses apart from buses 4, 
8, 9, and 12. Loads are fixed to 80% of their peak value, and 
reactive power marginal costs are set to c„ = 0.0125(];/kVar& 
h for all n. 
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Fig. 3. IEEE 13-bus feeder iMl . 
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Pig. 4. Voltages for buses 2, 12, the IEEE 13-bus feeder. 


Three reactive power control rules are tested; (i) the proxi¬ 
mal gradient descent (PGD) of (l24li : (ii) the proximal gradient 
descent with diagonal scaling (DPGD) in (l28] i: and the accel¬ 
erated proximal gradient descent (APGD) of (|29] |. The squared 
voltage magnitudes obtained at three buses are illustrated in 
Fig. a For all control rules, the step size is conservatively 
set to /i = 0.1/Aniax(X). Figure a demonstrates that APGD 
has an obvious four-fold speedup advantage over PGD and 
DPGD. Diagonal scaling does not exhibit any convergence 
rate advantage over PGD. The latter could be explained by 
the fact that the diagonal entries of X have similar values; 
hence, matrices X and D^/^XD^/^ with D = diag(diag(X)) 
have similar condition numbers. On a different note and as 
expected, the minimizer of (l22l) does not guarantee small 
voltage magnitude deviations across all buses: Buses 12-13 
exhibit small deviations, yet bus 2 converges to deviating from 
nominal by 1.03. 

The second experiment involves the SCE 47-bus grid with 
k(X) = 16,470 ifT^ . The capacitor located on bus 32 is 
ignored. Loads are set to 80% of their peak value with a power 
factor of 0.8. Five PV generators generating 60% of their 
capacity are located on buses 13, 17, 19, 23, 24. Distributed 
PVs with 50% penetration level are further installed on buses 
11, 12, 14, 22, 25, 33, 38, 39, 41. The relative cost value error 




Pig. 6. Voltages for buses 16, the SCE 47-bus feeder. 


attained by (l22l i using the optimal p = l/Aniax(X) is depicted 
in Fig. |5] Apparently, the novel scheme converges at least six 
to ten times faster than its alternatives. 

We further tested the accuracy of the linearized model over 
the full AC model calculated using the forward-backward 
algorithm ll22ll . Voltage magnitudes obtained from PGD and 
APGD with /i = 0.1/Aniax(X) are presented in Fig. |6l The 
curves suggest that the linearized model of (O offers a good 
approximation. 

To evaluate the control schemes on more realistic condi¬ 
tions, the third experiment uses the IEEE 123-bus grid and it 
also entails a topology reconfiguration; After twenty control 
periods (algorithm iterations), the switch between buses 97 
and 197 opens, while the switch between buses 151 and 300 
closes. Renewable (PV) generation units are located on buses 
32, 51, 64, 76, 80, 93, and 114, with capacities 60, 60, 120, 80, 
30, 100, and 80 kVA, respectively. The condition number for 
this feeder is «;(X) = 20,677. Figure |7] presents the squared 
voltage magnitudes over three representative buses and for 
fj, = 0.1/A„iax(X). Three observations are in order. First, note 
that the convergence guarantees for all three control schemes 
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Iterations 

Fig. 7. Voltages for buses 9, 105, and 149, for the IEEE 123-bus feeder. 



hold for any feasible initialization point. Theoretically, a 
topology change could potentially delay APGD since the 
parameters fit in (1291 ) are time-increasing. According to the 
curves, the topology change does not affect significantly the 
convergence rate of any of the algorithms. Second, the APGD 
scheme exhibits superior convergence properties over the PGD 
rule. Third, compared to the experiments on the IEEE 13-bus 
and the SCE 47-bus feeders, the longer convergence period 
can be attributed to the larger size of the feeder. 

B. Multi-Phase Grids 

To evaluate the findings of Section IIVI the control rule 
described by the iterations in (l3^ is tested on the multi¬ 
phase IEEE 13-bus system. Loads and PV penetration are 
selected as in the single-phase experiment presented earlier. 
The feeder voltage magnitude is fixed to 1, while the step size 
is set to n = A„ii„(A-i/2U^(X-f X^)UA-i/2). Squared 
voltage magnitude profiles obtained from (l38]) are plotted 
in Pig. [S] Voltages are calculated using both the linearized 
model and the full AC model at every iteration. Verifying the 
findings of ini, the curves indicate that the approximation 
is quite precise. The control rule converges within 40 itera¬ 
tions. Tests conducted with smaller step sizes exhibit slower 
convergence, whereas the scheme diverges for a step size 
M > 3.1A„,i„ (A”i/2U^(X -P XT)UA”1/2). 

VI. Conclusions 

To derive communication-free solutions, voltage regula¬ 
tion goals were lowered to penalizing large voltage mag¬ 
nitude deviations. Por balanced distribution grids, localized 
(un)constrained schemes were analyzed, while a control 
scheme based on proximal gradient descent was developed. 
Its computational complexity is comparable to that of the 
control rule suggested by the IEEE 1547.8 standard llT9l . yet 
it enjoys precise convergence guarantees. Adding memory to 
the control rule further yielded a significantly faster voltage 
regulation scheme. Por unbalanced distribution grids, a linear 
approximation model revealed a counter-intuitive inter-phase 


coupling. Based on the properties of the involved matrices, 
the developed reactive power control rule converges to a well- 
defined equilibrium point. Determining whether the latter point 
corresponds to the minimize!' of a particular optimization 
problem is a challenging open question. Numerical tests on 
benchmark feeders indicated the superiority of the acceler¬ 
ated scheme over diagonal scaling, and the resiliency of all 
novel schemes against topology reconfigurations. Incorporat¬ 
ing active power curtailment and performing local per-bus 
processing across phases constitute interesting future research 
directions. 


Appendix 

Proof of Proposition \2} Collecting 


Vtt — V = 2 Re 
(fTOl i. as 


bdiag({Z„})S 


for all n yields 


, or via the permutations of 


T(v^-v) = 2Re bdiag({Z„})TS 


(41) 


Pocus first on the LHS of (l4Tl l. Observe that voltage drops 
between adjacent buses can be alternatively expressed as 
V 7 r ,0 — V 0 = Av 0 + uoao for f S {a,b,c}. Stacking the 
latter equations across all phases yields 


v,r - V = (I 3 C) A)’V + W 0 I 3 ® ag. (42) 

Regarding the RHS of (HTt . recall that flows are decoupled 
across phases. The grid topology and (ISal i dictate that = 
A^S^, or equivalently, for all f. Stacking flows 

across all phases yields 

S = -(l3®F^)s. (43) 

Plugging (|42]|-(|4^ into (ED, and solving for v results in 

V = - 2 Re [(I 3 O A)~iT-i bdiag({Z„})T(l 3 O F^)s] 

- i;o(l3 0 A)~^(l3 (g) ao). (44) 


Using the facts A ^ = —F, T ^ = T^, Fag = Itv, and 
properties of the Kronecker product, (l44l i becomes 


V =2 Re 


(I3 0 F)T^ bdiag({Z„})T(l 3 0 F^)s 


■VqI 


3iV- 
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Substituting s = T^s and v = T'^v proves the claim. ■ 
Proof of Proposition |5} The claim is an application of 
[7, Prop. 2.1]. The Lipschitz constant of V/i(q) is A^ax(^)- 
Observe also that ||g ‘||2 = (v* — t;ol)^XX^(v* — uq!) < 
ALx(X)||v* -uollli, or 


|v - U 0 III 2 > 


1 


Amax(X) 


|g‘l|2. 


Note further that (v* —uol)^g* = (v* — wol)^X(v‘—riol) > 
Amin (X) 11V* — Wo 1II 2 ’ implying that (v* — wq!) is a descent 
direction for /i(q*) - although not the steepest one. ■ 

Proof for Proposition [j]- As in Prop. [2 showing Prop. |5] 
relies on an application of [7, Prop. 2.1]. In this case, the Lip¬ 
schitz constant for V/i(q) is Amax(X^X), and it also holds 
that ||vt - fiollb > Am[n^(X^X)||V/i(q*)|| 2 . The critical 
point in the multi-phase grid case though is that the positive¬ 
definiteness of Xj: guarantees that (v* — wqI) is a descent 
direction, since (v* — uol)^V/i(q*) = (v* — wol)^X^(v*— 

UqI) = (v‘ - Uol)^X^(v* - Wol) < Amin(Xa;)||v‘ - Wol|||. 


Proof of Proposition |6} By definition, it holds that 

\\I-pX\\l = Amax ((I - MX)^(I - mX)) = Amax (I " /^X^) 

where X^ := X+X^ —pXX^. Guaranteeing ||I —p,X ||2 < 1 
is equivalent to satisfying Amax (I — ftX^) < 1, or, simply 
1 — /iAmin (X^) < 1. Granted that p. > 0, the latter is 
equivalent to ensuring X^ to be a positive-definite matrix, 
i.e., X-bX^ /rXX^. It can be easily verified that pre/post- 
multiplying the aforementioned linear matrix inequality by 
yields the condition imposed on p by 
Proposition |6] ■ 
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